UCM/CSIC-95-07 
July 1995 



Real Space Renormalization Group 
Methods and Quantum Groups* 

Miguel A. Martm-Delgadof and German SierraJ 

\Departamento de Fisica Tedrica I 
Universidad Complutense. 28040-Madrid, Spain 
\Instituto de Matemdticas y Fisica Fundamental. C.S.LC. 
Serrano 123, 28006-Madrid, Spain 



Abstract 

We apply real-space RG methods to study two quantum group invariant Hamiltonians, that of 
the XXZ model and the Ising model in a transverse field defined in an open chain with appropiate 
boundary terms. The quantum group symmetry is preserved under the RG transformation except 
for the appearence of a quantum group anomalous term which vanishes in the classical case. We 
obtain correctly the line of critical XXZ models. In the ITF model the RG-flow coincides with the 
tensor product decomposition of cyclic irreps. of SU q (2) with q A = 1. 
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Real Space Renormalization Group methods, as applied to quantum many-body Hamiltonians, origi- 
nated from the successful study of the Kondo problem by Wilson Q. Later on people working in field 
theory and condensed matter generalized it to other problems by using the Kadanoff 's concept of block 
P], H. The Block method (BRG) has the advantage of being conceptually and technically simple, 
but it lacks of numerical accuracy or may even produce wrong results. For this reason the analytical 
BRG methods were largely abandoned in the 80's in favor of numerical methods such as the Quantum 
Monte Carlo approaches. In the last few years there has been new developments in the numerical RG 
methods motivated by a better understanding of the errors introduced by the splitting of the lattice 
into disconnected blocks. A first step was put forward in Q where a combination of different periodic 
boundary conditions applied to every block lead to the correct energy levels of a simple tight-binding 
model. This method however has not been generalized to models describing interactions. A further 
step in this direction was undertaken by White in || were a Density Matrix algorithm (DMRG) is 
developed. The main idea is to take into account the connection of every block with the rest of the 
system when choosing the states which survive the truncation procedure. The standard prescription 
is tho choose the lowest energy states of the block Hamiltonian. Instead, in the DMRG method one 
replaces the block Hamiltonian by a block density matrix and chooses the eigenstates of this matrix 
with the highest eigenvalues. The density matrix is constructed out of the ground state of a superblock 
which contains the desired block. 

In this letter we want to propose another RG method which uses the concept of quantum groups. 
This mathematical notion emerged in the study of integrable systems and it has been applied to 
conformal field theory, invariants of knots and manifolds, etc. ||, Q. The new application of quantum 
groups that we envisage has been partially motivated by the aforementioned work of White, Noack 
and collaborators an it is probably related to it. This relation is suggested by the fact that quantum 
groups describe symmetries in the presence of non-trivial boundary conditions. The typical example to 
understand this property of quantum groups is given by the ID Heisenberg-Ising model with anisotropic 
parameter A. The isotropic model A = ±1 is invariant under the rotation group SU(2), but as long 
as |A| 7^ 1 this symmetry is broken down to the rotation group U(l) around the z-axis. One can 
"restore" this full rotation symmetry by adding appropiate boundary operators to the Hamiltonian 
of the open chain. The classical group SU(2) becomes then the quantum group SU q (2), where the 
quantum parameter is related to the anisotropy by A = 2±| — ||, ||. The "restoration" of a classical 
symmetry into a (/-symmetry is achieved at the price of deforming the algebra and the corresponding 
addition rule of angular momentum. The (/-sum rule, which is called the comultiplication, becomes non 
local and violates parity. The total raising (lowering) operators acting on the whole chain are a sum 
of the raising ( lowering) operators acting at every single site times a non-local term involving all the 
remaining sites which appear in an asymmetric way: sites located to the left or to the right at a given 
site contribute differently. 

These features of (/-groups made them specially well-suited to implement a RG method which takes 
into account the correlation between neighboring blocks. Let us next show how this can be done 
explicitly in two examples in ID: the Heisenberg-Ising model and the Ising model in a transverse field 
(ITF). 



Heisenberg-Ising Model (XXZ Model). The open spin chain Hamiltonian is defined as: 

N-l 

H N (q,J)=J2hj,j+i(q,J) (!) 
h jJ+1 (q, J) = J -[o x 3 <?* +1 + ay j+1 + Z±fLa*aj +1 - ^f^ip*, ~ (2) 



1 



where ffj are standard Pauli matrices acting at the j site of the chain. For the time being q is an 
arbitrary complex parameter. Observe that the successive terms in cr| — <?f +1 in (||) when added into 
the total H only gives boundary operators proportional to a\ — a z N plus the standard bulk Hamiltonian. 
Let us now introduce the g-group generator of SU q {2) acting in the spin chain of ./V sites |J. 

1 N 

N 

S ± = ^2 g-a W + - +0 ?-i)o^ g 3(^+i+-+^) (4) 

3=1 

which satisfy the quantum group algebra, 

[S + ,S~] = q2S ~ Q ™ (5) 
q-q 1 

The important fact is that not only the whole Hamiltonian (||) but also the site-site Hamiltonian (|2|) 
commutes with the generators (^)-((5|) of SU g (2): 

[h jd+1 ,S*] = [h j , j+1 ,S ± ]=Q Vj (6) 

Hence the eigenstates of Hn(Qi J) can be classified according to the representations of SU q {2) (see || 
for details.) 

To construct a real space RG for the Hamiltonian (|l]) we shall choose blocks of 3 sites. This is 
important in order to get a renormalized Hamiltonian of the same form as the original one. The block 
Hamiltonian involving the first 3 sites is simply, 

H B = h 12 + h 23 (7) 

Now we can apply g-group representation theory to diagonalize Hb- There are 3 energy levels corre- 
sponding to the g-tensor product decomposition 1/2 ® 1/2 ® 1/2 = 1/2 © 1/2 © 3/2. 
For q real and positive (i.e. A > 1 ) or q a phase (i.e.|A| < 1) the lowest energy level of (|7|) is doubly 
degenerated and corresponds to one of the spin 1/2 irreps, which reads explicitly 



1. 1 



2' y / 2 (q + q^ + 



T (~q- 1/2 \ ITT) + (q 1/2 + q~ 1/2 )\ TIT) - q 1/2 \ TU>) (8) 



1 , 



„ (q 1/2 \ UT) " (q 1/2 + q- 1/2 )\ IU) + q 1/2 \ TIT)) (9) 



2' y /2{q + q-i + 

whose energy is es = — i(q + q~ l + 2). 

If we take the q — > + limit in (|l]) and (Q) we obtain an Ising model Hamiltonian (A — > oo) with 
a unique ground state given by the Neel state | ITIT • • •) (notice that this uniqueness is due to the 
boundary term — (o-f — a z N )). On the other hand, for a block of 3 sites there are 4 states of lowest 
energy (| TIT); I ITT)> I ITI) an d| IIT)); while for a block of 4 sites there is again only one ground state 
given by | HIT). This means that choosing an odd number of blocks is not appropiate to study the 
Ising limit of (§). To do so one should choose an even number of sites, but this will not be pursued 
here. Hence we shall concentrate on q being a phase. 

The renormalization prescription consits in choosing the states as the spin up | T)' and down 

| I)' states associated to the whole block as if it were a single site. Using the standard methods of BRG 
we obtain the following RG-transformation laws for the spin operators S{ acting at the sites i = 1 and 
3: 



2 



(Sf) RG = ^q) S' x a = l,3 (10) 

(s!) RG = aq)s' y * = i,3 (ii) 

(Sf)BG = e(?)fi ! '* + 7 W ( g )l / . = 1,3 (12) 
where £(q) is a renormalization factor which depends upon q as: 

«9) = ? + g " 1 1 + \ (13) 

and 

= = ^) = 4( / + ~?i + 1) ( 14 ) 

The multiplicative renormalization factor £(q) is common to all the spin operators <Sj as a consequence 
of the full symmetry group SU q {2). The "quantized" feature of SU q (2) is reflected in the "quantum 
group anomaly" term in (|l2])-([l^), which indeed shows the deviation from the classical case (q = 1). 
Eqs. ([lO|)-( |T2]) a re quite different from the standard BRG analoge for the Heisenberg-Ising model done 
in reference PQ l, where the RG-equations for S x and S y differ from those of S z (i.e. £ x = £ y ^ £ z .) 
Using eqs. (|To|) -(12) we can get the renormalized block-block Hamiltonian h 3 k,3k+i- Putting all terms 



together we arrive at the following effective Hamiltonian H' which acts on the chain having N/3 sites, 



H' = H N/3 (q', f) + je B (q, J) + (j - 1) e BB (q, J) (15) 



where 



Q = Q (16) 
J' = i\q)J (17) 

e B (q,J) = -^{2 + q + q- 1 ) (18) 

e BB(q,J) - J 16{q + g -l + 1) 2 (I 9 ) 

The e B contribution to the energy comes from the block part Hb while bbb is a novel contribution 
coming from the quantum group anomaly. The remarkable feature of eqs. (15) and (|lfj) is that the 



coupling constant q, or alternatively A does not flow under the RG-transformation, while J^ m \ which 
is the value of J after m RG-steps, goes to zero in the limit where m —* 00, which in turn implies that 
the theory is massless. Hence, eqs. (|l6|) predicts correctly a line of critical models in the range | A| < 1. 
These models are described by Conformal Field Theories (CFT) with central extension c less than one. 
If we write the quantum parameter as q = e J7r /^ +1 ^ then c = 1 — jjjnfjr^n I P- 11 - The boundary terms in 
(|l])-(Q) are responsible for this fact. 

A non-trivial check of the validity of our RG-method can be given in the case where q = e l7r//3 for 
which c = (percolation limit). The ground state energy can be computed exactly from the constant 
terms ([TJJ) (assume that N = 3 m and perform m-RG steps), and is given by 

Eo (N,q = e i ^ 3 ) = ~N + ^ (20) 
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This equation coincides with the exact result obtained in |11] using Bethe ansatz. This means that at 
least in this particular case the truncation to the states | ± A) in dH)-(^) of our gRG-method does not 
imply any kind of approximation. This is consistent with the fact that in the CFT with c = there is 
a unique state, namely, the ground state. What the gRG-method does is to pick up that piece of the 
ground state which projects into a given block. 

It can be shown that a consistent representation theory of quantum groups at root of unity (i.e., 
q^ +1 = —1) requires the use of truncated tensor products of g-group irreps. In the case of q = e ln ^ 
this truncation implies: 

(1/2 ® 1/2 ® l/2) ?=eW3 = 1/2 (21) 
which is precisely the truncation performed when restricting ourselves to the states (@)-@. In ref- 



erences [12] the representation theory of g-groups was put in one-to-one correspondence with that of 
Rational Conformal Field Theories (RCFT). There it was observed that the truncation inherent in the 
construction of the RCFT's has a parallel in the truncation of the representation theory of g-groups 
with q a root of unity. The result we have obtained in this letter suggests that g-group truncations can 
be carried over a RG analysis of g-group invariant chains. In other words, using g-groups we can safely 
truncate states in the block RG method. We may summarize this discussion squematically by, 

gRG-truncation <-> RCFT (22) 



Ising model in a transverse field (ITF). This simple model has been widely used to test the validity 
of BRG methods Q, Q. The Hamiltonian of an open chain is given by H = J2f=i~ where 

h jd+1 = -(Ja*o-* +l + pa] + pV| +1 ) (23) 

The standard choice is p = p' = T/2, in which case ( p3|) has 4 different eigenvalues. The BRG method 
with a block with two sites chooses just the 2 lowest ones. However if (p,p') = (T, 0) (or (0,T)) the 
Hamiltonian ( |23"| ) has two doubly degenerate eigenvalues ±e_B (cb = V J 2 + T 2 ). This choice is not 
parity invariant but it implements the self-duality property of the ITF model, yielding the exact value 
of the critical fixed point of the ITF which appears at (T/J) c = 1 filjfl . This degeneracy of the spectrum 
of ( |23"| ) has a q-group origin. The relevant quantum group is again SU q {2) with q — 1. However the 
representations involved are not a (/-deformation of the spin 1/2 irrep. as in the previous example, but 
rather a new class of irreps. which only exist when q is a root of unity. They are called cyclic irreps. 
and neither are highest weight nor lowest weight representations as the more familiar regular irreps. If 
we call E, F and K the generators of SU q (2), which correspond essentially to S + , S~ and q 2S * in the 
notation of the previous example, then a cyclic irrep. acting at a single site of the chain is given by: 

Ej = oaf, Fj = ba y j, Kj = \a] (24) 

where a = ^VA 2 — 1, b = — |vl — A~ 2 . The parameter A is the label of the cyclic irrep. Using (|24|) 
and the addition rule of SU q {2) we can get the representation of E, F and K acting on the whole chain: 



N 



3=1 
N 



(25) 
(26) 



4 



N 

K = X N l[a z j (27) 

3=1 



Now it is a simple exercise to check that these operators commute with (|23|), 

[h jtj+1 , E] = [h jtj+ i, F] = [h jij+1 , K] = 0, Vj (28) 



assuming that we choose 



A = T/J (29) 



The last of the equalities in (^) expresses the well-known ^-symmetry of the ITF- model which allows 
one to split the spectrum of the Hamiltonian into an even and odd subsectors. The other two symmetries 
are new and explain the degeneracy of the spectrum of hjj+i. By all means the whole Hamiltonian 
H = hjj + i is also invariant under (25). Notice that H differs from the standard ITF simply in 



a term at one of the ends of the chain. This is the same mechanism as for the XXZ Hamiltonian: 
one needs properly chosen operators at the boundary in order to achieve quantum group invariance. 
Similarly as for the XXZ model the RG-analysis of the ITF becomes a problem in representation of 
quantum groups: blocking is equivalent to tensoring representations. What is the tensor product of 
cyclic irreps.? Here it is important to realize that all cyclic irreps. of SU q {2) have dimension 2, what 
distinguishes them is the value of A. The tensor product decomposition of two cyclic irrep. Ai and A2 
is given by: 

[Ax] ® [A 2 ] = 2 [A1A2] (30) 

where the 2 means that A1A2 appears twice in the tensor product. If we perform a blocking of two sites 
we will get two cyclic irreps. corresponding to A 2 . Then we expect from (/-group representation theory 
that the new effective Hamiltonian hj j +1 will have the same form as ( p3|) but with new renormalized 
coupling constants J' and V satisfying: 

A' = j, = (j) 2 = A 2 (31) 



This is indeed the result obtained in [|T3|]. We arrive therefore at the conclusion that the RG-flow of 
the ITF Hamiltonian $2^ ) is equivalent to the tensor product decomposition of cyclic irreps of SU q {2). 
This g-group interpretation of the RG-flow is independent of the size of the blocks: for a n-site block 
the RG-flow would be A — > X n . The fixed point A = 1 of ( |3l| ) describes the critical regime of the ITF 
Hamiltonian and it corresponds to a singular point in the manifold of cyclic irreps. |lq] , |i~5fl . At A = 1 
the operators ( [25] ) are still symmetries of the Hamiltonian (a, b taking any non-zero value) and they 
recall the Jordan- Wigner map between Pauli matrices and ld-lattice fermions. 



Cyclic irreps. were used in [16] to derive the Boltzmann weights of the iJ/v-chiral Potts model [17]. 
In [|r| the labels of the cyclic irreps. have the meaning of rapidities rather than coupling constants as 
in our realization. 

The Z2 CP-model is nothing but the ITF model. We may wonder whether the general ^jy-chiral 
Potts model admits a g-group RG treatment along the lines of this work. This problem will be considered 
elsewhere. A model that admits a qRG analysis is the XY model with a magnetic field h (XYh). The 



results will be presented in [18]. It suffices to say here that the q-group underlying the model is SU q {2) 
with q 4 = 1 and the representation used are the so called nilpotent irreps. [j?]], which are also described 
by a parameter A analoge to that in ( f24|) and related to the magnetic field h. The XYh model is 
equivalent to a free fermion with chemical potential. The results we are obtaining can be translated 
into a g-group symmetry between fermions, either free as in the XY or ITF models or interacting as 
the XXZ model. Another interesting model of interacting fermions is the Hubbard model, which has 
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been studied using RG-methods in |19[ . The integrability of the ID Hubbard model [|2G|] suggests that 
it might be studied using our qRG techniques. 

All the Hamiltonians analysed in this letter are one-dimensional, so they are quantum groups as we 
know them. Despite the fact that the Yang-Baxter equation (the precursor of g-groups) has a higher 
dimensional analogue called the Zamolodchikov or tetrahedron equation [21], the corresponding high 
dimensional analogue of quantum groups is not known. This fact represents a barrier to a qRG analysis 
of Hamiltonians defined in dimensions higher than one. 

Another possibility, which is suggested by our results, would be to define quantum groups as those which 
contain symmetries which are anomalous under RG transformations. This definition is independent of 
the space dimensionality. The quantum anomalous term in equation ( |l2] ) gives a discrete realization 
of this idea. A continuum analogue of this anomaly is given by the Feigin-Fuchs current, which has 
an anomalous operator product expansion with the energy- momentum tensor [22]. At this point it 
may be worth to recalling the continous version of quantum groups in CFT of reference |23|| , which 
uses the Feigin-Fuchs or free field realization of the latter. Putting all these arguments together, 
we arrive at the conclusion that quantum groups are indeed defined by symmetries anomalous under 
RG transformations. This point of view about quantum groups may set up the pathway to new 
developments in the field. 
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